Geometrical barriers and the growth of flux domes in thin ideal superconducting disks 
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When an ideal (no bulk pinning) flat type-II superconducting disk is subjected to a perpendicular 
magnetic field Ha, the first vortex nucleates at the rim when Ha = Hq, the threshold field, and 
moves quickly to the center of the disk. As Ha increases above Ho, additional vortices join the 
others, and together they produce a domelike field distribution of radius b. In this paper I present 
analytic solutions for the resulting magnetic-field and sheet-current-density distributions. I show 
how these distributions vary as b increases with Ha, and I calculate the corresponding field-increasing 
magnetization. 

PACS numbers: 74.25.-q,74.78.Bz,74.25.0p 
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I. INTRODUCTION 

The penetration of magnetic flux into noneUipsoidal 
superconducting samples is complex. The initial en- 
try of magnetic flux in the form of multiply quantized 
flux tubes in type-I superconductors or singly quan- 
tized vortices in type-II superconductors is impeded both 
by a barrier of geometric origin and by bulk pinning. 
The effects of the geometrical barrier are most pro- 
nounced in type-I superconductorsii^i^i^ and in type- 
II superconductors in which the bulk pinning is unusu- 
ally low.— iii^ Detailed studies of geometrical-barrier ef- 
fects in strips and long slabs, in which the magnetic- 
field distribution is two-dimensional and conformal- 
mapping techniques can be applied, have been car- 
ried out for both type-I superconductor o^'^° and type-II 
superconductors^iaiia2i3a4i5^a7 

Various authors have carried out numerical calcula- 
tions of the magnetic-field and current distributions in 
flat disks of various radii and thicknesses, including the 
possibility of the formation of a cylindrically symmetric 
flux dome at the center When the flux 

distribution in the disk depends only upon the radial co- 
ordinate p, the problem is two-dimensional, in that the 
magnetic field can be assumed to depend only upon p 
and the axial coordinate z. However, conformal-mapping 
techniques can no longer be applied in this geometry. 

In this paper I present analytic solutions for the 
magnetic-field and sheet-current-density distributions for 
the geometrical-barrier problem in thin superconducting 
disks. As in most of the previous theoretical studies, 
I have ignored the spatial structure on the scale of the 
spacing between fiux tubes or the intervortex spacing and 
have considered the magnetic flux density averaged over 
this fine scale, such that the solutions have cylindrical 
symmetry. Although I describe magnetic-flux penetra- 
tion in terms of the entry of vortices into type-II super- 
conducting disks, the theory could be applied with only 
minor modifications to the entry of fiux tubes into type-I 
superconducting disks. 

The paper is organized as follows. In Sec.|TTl I consider 
the initial entry of vortices into a disk and the formation 



of a fiux dome whose radius h is much smaller than the 
disk radius a. In Sec. IIIII I theoretically examine the 
subsequent entry of magnetic fiux as the dome radius in- 
creases and h approaches a. In Sec. IIVI I present results 
for the field-increasing magnetization, and in Sec. |Vl I 
briefly summarize the main findings of this paper. Four 
Appendices contain numerous mathematical details un- 
derlying the theory. Finally, Sec. lVII contains some of my 
recollections of John Bardeen, who is being honored in 
this issue. 



II. INITIAL GROWTH OF SMALL FLUX 
DOMES 

Consider an ideal (pinning-free) thin type-II supercon- 
ducting disk of thickness d and radius a centered on the 
z axis in the xy plane. I assume that either a ^ d 
if d/2 > A, where A is the London penetration depth, 
or a > A if d/2 < A, where A is the Pearl length)^ 
A = 2X^/d. In the Meissner (vortex-free) state, when a 
magnetic field Ha is applied along the z direction, the 
magnetic field in the plane z = is, as shown in Refs.[23l. 
m, and [13 and Appendix El 



Haz{p,0) = 0, p < a, 

r 2 

p > a. 
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When p ^ a, Haz{p,0) is dominated by the term with 
the inverse square-root singularity. 

The corresponding sheet-current density in the disk 
(p < a) is, from Eq. (|Al6p . 
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Using the reasoning of Ref. |23|, I make the simphfying 
assumption that the first vortex enters when the applied 
field reaches the value Hq, at which the demagnetization- 
enhanced field at the edge, estimated as 



(5) 



where 5 is the larger of d/2 or A, becomes equal to Hg, 
where Hg = Hd, the lower critical field, if there is no 
Bean-Livingston barrier, or Hs ~ He, the bulk thermo- 
dynamic field, if the edge is without defects and thermal 
activation is negligible. Thus 



(6) 



However, to determine precisely when a vortex en- 
ters or exits at the edge is a diflacult prob- 
lem, because these processes are sensitive to de- 
tails such as the shape and perfection of the 
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When Ha is just above Hq, vortices enter a pin- free 
disk and collect in a flux dome at the center, driven by 
the Lorentz force, 



Close to the center of the disk. 



-AHa4>ap/Tra. 



(7) 
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As shown in Fig.[Tl the vortices in the flux dome gener- 
ate a magnetic field z) = pHyp{p, z) + zHyz{p, z), 
which is calculated in Appendix [B] The return flux 
around the rim of the disk partially cancels the fleld 
Ha,edge- The vortex-geucratcd local field at the outer 
edge can be estimated as [see Eq. (jB2|) ] 



Hv,edge ~ Hyz{a + 5, 0) 



(9) 



The condition for the next vortex to enter the strip is 
therefore that 



H 



edge 



^H„ 
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v,edge 



^H. 



or 



Ha — Ho — 



2'Kpf)0? ' 



(10) 
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such that the total magnetic flux in the fiux dome of 
a superconducting disk initially grows linearly with the 
difference Ha — Hq. 

Next, let us examine how the radius h of the dome 
initially varies as a function of Ha — Hq. The dome- 
generated magnetic field near the center of the disk looks 
much like that produced by a distribution of magnetic 




FIG. 1: Contour plot of the stream function ipv{p,z), where 
the contours correspond to lines of magnetic field emanating 
from a flux dome carrying flux $ up through the disk at the 
origin. The flux <I> returns around the perimeter of the disk 
at p = a = 1. 



monopoles. When 6 <C a, the dome's magnetic-field dis- 
tribution can be calculated to good approximation by 
treating the problem as if the disk were of infinite radius. 
As shown in Appendix [Cl the magnetic flux $ penetrat- 
ing up through the strip spreads out within p < b with 
the distribution Hsz{p,Q) = Hd^Jl - p^ /h"^ . The corre- 
sponding dome-generated magnetic field tangent to the 
top surface within p < b for a small dome is 



Hspip,0- 



46 ■ 



(12) 



The field on the bottom surface has the same behavior 
but of opposite sign, such that the dome-generated sheet- 
current density behaves as 

T^HdP 



Ks4>{p) 



2b 



(13) 



However, the net sheet-current density K^{x) inside the 
dome must vanish in order that there be zero net Lorentz 
force on any vortex within the dome. Setting K^{p) = 
Ka(t>{p) + Ks(t){p) = yields for a small dome 



H, 



;Ha- 



(14) 



Outside the dome, assuming b < p a, the results of 
Appendix [C] can be used to show that 



i^.(p) = -i^.[f(f 



1 b 



— sm — 



62 



p 



(15) 



Since K^{p) < 0, the radial Lorentz force on any vortex 
outside the dome obeys Fp{p) — K^{p)(j)o < 0, such that 
all vortices are confined within the dome. 
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FIG. 2: Dome radius vs increasing applied field, expressed 
as h/a vs Ha/Ho. The dashed curve shows the small-dome 
approximation of Eq. (|18p , and the solid curve shows the the- 
oretical results for domes of arbitrary radius, obtained from 
Eq. dMj. 



The magnetic flux contained within a small dome is 

$ - 27rAio / H,,{p)pdp - 2T:poHdby3. (16) 
Jo 

Combining this result with Eqs. pTjl and leads to 
the condition 



Ha — Hq 



0.270(-) , 



(17) 



III. ARBITRARY FLUX DOME RADIUS 

As the applied field Ha increases and additional vor- 
tices move to the center, the radius h of the dome ex- 
pands, and the small-dome approximations used in the 
previous section are no longer valid, because they as- 
sume that b <^ a. In particular, the small-& dome 
shape, Hsz{p,0) = Hdy/l — p'^/b'^, is no longer vaHd. 
Appendix |D] contains the details of the magnetic-field 
and sheet-current distributions resulting from an array 
of vortices carrying magnetic flux $ up through an an- 
nulus of radius pc < a. These results can be used to 
obtain the Green's functions that generate the magnetic- 
field and sheet-current-density distributions produced in 
response to a large flux dome characterized by an arbi- 
trary field distribution Hj_z{p',0) confined to the region 
p' < b < a. We need only to replace pc by p' and <& by 
2t: fioHdzip' ,0)p' dp' and then integrate over p' from to 
b. The resulting expression for the z component of the 
dome-generated magnetic field for p > a is 



Hdz{p,0) = -- [ 

71" Jo 



Hdzip',0)p'Va^-p'^ 



This expression has the property that 



dp'. (21) 



In I 

J a 



Hdzip, 0)pdp = -2tt / Hdz{p, 0)pdp, (22) 



or 



which is equivalent to the statement that the return flux 
through the space p > a is the negative of the magnetic 
flux up through the dome p < b. Similarly, the resulting 
expression for the dome-generated sheet-current density 



3tt^\^/3 /'Ha- Hq\^/3 

Ha 



1.547 



Ha-H o\U3 
Ha 



(18) 

when b/a ^ 1. Numerical calculations of the dome ra- 
dius carried out as in Ref. [2^ are in agreement with these 
results. The dashed curve in Fig. [5] shows a plot of b/a 
given by the small-dome approximation of Eq. (jl8p for 
an increasing applied field Ha- 

If the applied field Ha is reduced after a small dome 
has appeared at the center, the magnetic flux $ in the 
dome remains constant because the edge conditions no 
longer permit the entry of vortices, and since K^{p) < 
for & < p < a, any vortex that strays out of the dome 
is forced back in. Thus as Ha decreases, the dome ra- 
dius b expands and the field Hd at the center decreases. 
According to Eqs. ((H]) and (fTB)) . 



b=( 



37r<I>a \i/3 



V 16fJ.oHa 



(19) 



and 



KM = 2 Hdz{p',Q)g{a,p',p)p'dp', (23) 
Jo 

where Q{a,p',p) is defined in Appendix ID] 

Because the disk is in a perpendicular applied field 
Ha^ which induces its own sheet-current density Ka,j,{p) 
given in Eq. ([3]) , the net sheet-current density is (p) = 
Kacjiip) + Kd<p{p)- For the vortices in the dome to be in 
equilibrium, the net Lorentz force on any vortex within 
the dome must obey Fp{p) = K^{p)(j)o = 0. The field dis- 
tribution Hdz [p, 0) inside the dome therefore must obey 
the integral equation 



Hdz{p',0)g{a,p',p)p'dp' 



2HaP 



(24) 



for all p < b. 

Another equivalent integral equation can be obtained 
by noting that the condition K^{p) = 2Hp{p, 0+) — for 
p < 6 is equivalent to the condition that d4>{p, 0+)/dp — 
or 



(20) 



0(p,o+) = Mp,o+) + Mp,o+) 

= 0d(O,O+) + 0a(O,O+), (25) 



4 



where 4'a{p,0+) = 4'a2{p,0+) is given in Eq. (jA6p . and 
4>dip,0+) can be obtained from (pc{p,0+) in the same 
way that Kii4,{p) was obtained from Kcc/,{p): 



Mp,0+)^ / Hd:,{p',0)T{a,p',p)p'dp'. 



(26) 



With the help of Eq. (jD16|) . Eq. (^5]) can be expressed 
as the following integral equation, from which Hdz [p, 0) 
can be determined: 



2Ha 



1 _ ^(a. p', p)p' dp' 
a i 

(a - V^^^), p<b. (27) 



Equation (|24p can be recovered by differentiating Eq. P?)) 
with respect to p. These equa tions are solved exactly 
by Hdz{p, 0) = {8/TT'^){b/a)Hay/l- p'^/b^ in the limit as 
b/a — > and by Hdz{p, 0) = -ffa when b = a. However, I 
have been unable to find analytic solutions for H^z {p, 0) 
for arbitrary values of b/a, unlike the situation for a flux 
dome in a superconducting strip, for which the integral 
equation can be inverted exactly for all dome sizes. 

On the other hand, numerical solutions of Eq. (l27l) 
reveal that the functional form of Hdz{p,0) is within a 
few percent of 



Hdz{p,i)) = Hd 



'i_fPjh^ 

l-p^/a^ 



(28) 



similar to the behavior of the dome shape in a 
strip,— except that here the value of Hd increases from 
{8 /tt'^) (b/a) Ha for b/a < 1 to iJ^ as b/a ^ 1. Thus, 
the value of Hd can be obtained to good accuracy by 
substituting Eq. (|28|) into Eq. ((27|) and evaluating it nu- 
merically at p ~ b. The results for Hd vs f3 — b/a for 
< /3 < 1 can be fit to 0.1% by 



Hd^HahiP), 



(29) 



where 



/i(/3) = /?(_ + 0.165/?^ - 0.028/3-* + 0.052/3''). (30) 

Recall that the growth of flux domes when the dome 
radius b obeys f3 = b/a 1 is determined by Eqs. (O and 
PH)) . For larger domes, the dome-generated local field at 
the outer edge can be obtained with good accuracy from 
Eqs. (HI]) and ^ as 



H. 



d,edge 



Hdz{a + S,0) 
2Hd p'^ 
26 Jo 




p'2) 



-dp' 



where 



/2(/3) = l-^^i^tan-i 



(31) 



(32) 



The condition for the next vortex to enter the strip is 
therefore that 



Combining Eqs 
for Ha > Hq 



Hedge ~ Ha,edge H~ Hd ^dge — Hg. (^3) 

29]), jH]), and ([33]), we obtain 



Ha 



1 



Ho 1 -<?(/?)' 



where (3 = b/a and 



5(/5) = /i(/5)/2(/3). 



(34) 



(35) 



A plot of the reduced dome radius (3 = b/a vs Ha/Ho, 
obtained from Eq. ([34]), is plotted as the solid curve in 
Fig. [2 

From Eqs. ((28]) and (I29|) we obtain the corresponding 
approximation to the total magnetic flux up through the 
disk into the dome, 



$ = 2^ / Hdz (p, 0)pdp = Hdh {P) = Hah W) h (/3) , (36) 



where j3 — b/a and 



HP) 



1- (1-/52) 



tanh ^/3 



P 



(37) 



Consider the situation when an increasing applied field 
reaches the value Hai, producing a dome of magnetic 
flux $1, radius 6i, and field Hdi at the center, where 
these quantities are related via Eqs. (|28]) . (|29]) . (|36]) . and 
([57]) . Suppose that the applied field Ha is now reduced 
below Hal- The magnetic flux $ in the dome remains 
constant at the value $i because the edge conditions no 
longer permit the entry of vortices, and since K^{p) < 
for b < p < a, any vortex that strays out of the dome 
is forced back in. Because both /i(/3) and fsiP) are in- 
creasing functions of /3, as Ha decreases, the dome radius 
6 expands and the field Hd at the center decreases. The 
new values of P and Hd can be obtained from Eq. ([M]) . 



IV. FIELD-INCREASING MAGNETIZATION 

The above results can be used to calculate the mag- 
netic moment per unit volume, Mz — mz/na'^d, which I 
shall refer to as the magnetization, where the magnetic 
moment is^ 



mz^TT I p K^{p)dp. 
Jo 



(38) 



Because K^p) = 2Hp{p,0+) = -2d(j){p,0+)/dp and 
(j){a, 0+) = 0, partial integration yields 



= 47r / p(l>{p,0+)dp, 

10 



(39) 
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where 0(p,O+) = (f>a{p,0+) + (j)d{p,0+) and the corre- 
sponding contributions to are Maz and M^^- 

Since 4>a{p,0+) = (/)a2(/0, 0+), given in Eq. (|A6p . we 
easily obtain the negative diamagnetic contribution aris- 
ing from the Meissner response to the appUed field, 

Maz = -Xo-ffa, Xo - 8a/37rd. (40) 

A positive contribution M^z due to vortices in the flux 
dome is obtained with the help of Eqs. and (|28p: 



f '' / sin 9 

Mdz = SxoHd / sin W 1 - —^n{0)d0, (41) 

Jo V sin 

where 9b = sin^^(6/a), 

ni9) = £;(sin6') 

_ 2 cos g[£;(cos 0) - sin^ 0K(cos ,^)] ^^^^ 
Jo TT cos (/)(1 — sin^ 6* cos^ 0) ' 

and K{k) and E{k) are complete elliptic integrals of the 
first and second kind of modulus k. Since 7i(0) = 1, 
when /3 = Va « 1, Mdz » Xo^Td/?' = XoHa{8/TT^)p^ = 
3xo{Ha — Hf)). In the hmit as b/a ^ 1, 9b ^ 7r/2, Hd 
Ha, the integral in Eq. (|4T|) yields 1/3, and M^z ^ XoHa- 
The total magnetic moment per unit volume Mz — 
Maz + Mdz, numerically calculated as above with the 
help of Eqs. dSQl), and is shown as the sohd 

curve in Fig. [31 Just above Hq, Mz « Xo(2^^a ^ Si/p), as 
shown by the dashed line. In the limit as Ha/Ho oo, 
Mdz — > XoHa, nearly canceling Maz, such that Mz 0. 
However, all the above analysis is valid only when Ha < 
Hc2, so that the superconductor remains in the mixed 
state. 

Suppose the increasing-field magnetization follows the 
solid curve in Fig. [3] in an increasing field Ha up to the 
value Hal, producing a dome containing magnetic flux 
$1. As described at the end of Sec. IIIIl if the applied 
field Ha is then reduced, the magnetic flux $ remains 
constant, but (3 — b/a increases and Hd decreases, and 
Eq. (15^ can be used to obtain both f3 and Hd as func- 
tions of Ha- In turn, these values can be substituted into 
Eq. ([iT|) to obtain Mz = Maz + Mdz as Ha decreases be- 
low Hal- This procedure has been used to calculate the 
resulting dot-dashed portion of the magnetization curve 
shown in Fig. [31 where Hai = 1.5Ho has been used as 
an example. This portion of the magnetization curve 
is reversible, provided that Ha changes slowly and Ha 
remains above the critical exit field H^^^, where b is suf- 
ficiently close to a that vortices escape from the sample. 

V. SUMMARY 

In this paper I have presented analytic calculations de- 
scribing the evolution of the flux dome of radius b pro- 
duced at the center of a thin ideal type-II superconduct- 
ing disk of radius a when an increasing perpendicular 




Ha/Ho 

FIG. 3; Magnetic moment per unit volume in units of 
XqHq vs increasing applied field Ha in units of Ho calculated 
from Eqs. (|40|l and (|41[) (solid). The short dashed line shows 
the behavior Mz — —xo{3Ho — 2Ha), which holds just above 
Ho when b/a <^ 1. The dot-dashed curve shows the reversible 
magnetization achieved after the applied magnetic field is first 
increased to 1.5-ffo and then decreased; the magnetic flux in 
the dome remains constant as Ha decreases and the dome 
radius b increases, thereby decreasing the area that can carry 
a screening current. 

magnetic field Ha is applied. This phenomenon, not ob- 
servable in samples with strong bulk pinning, is due to 
the geometrical barrier. The theory is exact in the lim- 
its b/a <C 1 and b/a — > 1. For intermediate values of 
b/a, this theory is approximate, but the results should 
be accurate to within one or two percent. 

I have found that b/a varies as {Ha — iJo)^^'^ when 
Ha is just above the threshold field Hq, but I also have 
shown how to calculate b/a with good accuracy for all 
values of Ha above Hq. In addition, I have shown how to 
obtain the magnetic field Hd at the center of the dome 
corresponding to each value of (3 = b/a. Both b and Hd 
should be accessible experimentally using scanning tech- 
niques with high field sensitivity and spatial resolution. 

I also have presented calculations of the magnetiza- 
tion (or more precisely, the magnetic moment per unit 
volume) of an ideal type-II superconducting disk as a 
function of Ha when the magnetization is dominated by 
the geometrical barrier. This quantity is also readily ob- 
servable experimentally. 



VI. JOHN BARDEEN 

Because this issue celebrates the 100th anniversary of 
John Bardeen's birth, Vladimir Kresin has asked me to 
add some personal recollections. When I was in my final 
year as a graduate student at the University of Illinois 
in Urbana, Leo Kadanoff, my first advisor, went off' on 
sabbatical leave and left me in the care of John Bardeen. 
I periodically reported my progress to Professor Bardeen 
and remember well that wonderful day when he said these 
precious words, "I think you have enough here for a the- 
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sis." I diligently went to work writing my Ph.D. thesis on 
anisotropy of the superconducting energy gap, and as I 
completed the chapters, I sent them to various scientific 
institutions around the world where John Bardeen was 
visiting, so that he could peruse what I had written. (By 
the way, my usually tolerant wife Judy reached the end 
of her patience, phoned me at my office, and instructed 
me to come home when I was working even on the 4th of 
July.) 

Although John had a quiet, modest manner, he did 
enjoy entertaining at his home at parties, which I expect 
were mostly organized by his wife Jane, an outgoing, gra- 
cious hostess. At such a party I recall a funny moment, 
which John also greatly enjoyed, when a fellow grad stu- 
dent Richard Craig ran around the back yard trying to 
catch a rabbit. 

John Bardeen was a soft-spoken man, and as students 
taking classes from him, we knew we had to sit in one of 
the front rows so that we could hear the pearls of wisdom 
he could pass on to us. I also remember when we grad 
students would sit through a long, obscure theory semi- 
nar given by a visitor, only to hear John Bardeen, during 
the question and answer period at the end, summarize 
in one clear sentence the speaker's main conclusions and 
why they were important. 

In addition to his brilliance as a scientist, John 
Bardeen had an unusual trait that it took me some time 
to get used to. In many conversations I have with others, 
our sentences are overlapping, with one person starting 
to speak before the other has finished. I suppose this 
is because most of us process input and output simulta- 
neously. (Of course, the disadvantage of this conversa- 
tional mode is that we too often speak before thinking.) 
I eventually learned that it was not this way with John 
Bardeen. When I said something to him or asked him 
a question, he took some time to think about his reply. 
While he was thinking, his input channel was turned off, 
and if I spoke during this interval, he could not hear me. 
I eventually learned to wait patiently for his reply, after 
which our conversation would continue. 

It was a privilege to have known John Bardeen, whose 
scientific contributions have had an enormous impact on 
both science and technology, and I appreciate this oppor- 
tunity to help celebrate his 100th birthday. 
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APPENDIX A: MEISSNER RESPONSE OF A 
DISK IN A PERPENDICULAR FIELD 

This Appendix contains a detailed description of the 
magnetic field in the vicinity of a superconducting disk 
of radius a centered at the origin when it is subjected 
to a perpendicular applied field Ha- In the absence of 
the disk, a uniform magnetic field in the z direction pro- 
duced by distant sources can be represented by a po- 
tential (paiiz) — —HaZ and a stream function ipaiip) = 
—HaP^/2. The magnetic field Hai = pHaip + zHaiz can 
be expressed as Haiz = —d(j)ai/dz — —{l/ p)d4'ail dp = 



Hn. and H, 



alp 



-d^alldp^{l/p)di^al{p)ldz = Q. 



The superconducting disk responds with induced 
screening currents, which generate a magnetic field that 
cancels the z component of the applied field inside the 
disk. The magnetic field resulting from the induced 
screening currents, Ha2 = pHa2p + zHa2z, is described 
by the potential 0a2(p, z) and stream function il'a2{p, z) 

via Ha2z = ~d(f>a2/dz -{1/ p)dlpa2/dp, and Ha2p = 



-d(j)a2/dp = {l/p)diJa2{p)/dz, whcrC V^(l)a2 



0. Here 



and in the rest of the Appendix, I use a formulation and 
method described by Lamb4S The solutions for z > are 



(I)a2{p,z) 
1pa2{p,z) 
Ha2z{p,z) 



fa(k)Jo{kp)e-'''dk, 

oo 

fa{k)Jo(kp)e-'''kdk, 
-^''kdk. 



Ha2p{p,z) = / fa{k)Ji{kp) 

Jo 

f (/fc) — 2i?Q d / sinfca 
TT dk\ k 



In the plane 2 = 0, 



2Hn 



<Aa2(p,0-l-) = -v^a^ - P^, P < a, 

TT 

= 0, p > a. 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 

(A6) 
(A7) 

(A8) 



sm 



P 

P>a, 



Ha2z{p,0) = -Ha, p<a, 

TT V 



Ha2pip,0+) = - 



p> a, 
2HaP 



— sm 



p- 



p<a, 



0, p> a. 



(A9) 
(AlO) 



(All) 
(A12) 

(A13) 



The overall field distribution, including both the ap- 
plied field and the Meissner response of the disk, is 
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Ha{p, z) = Hal + Ha2{p, z) = pHap{p, z) + zHaz{p, z) , 

which is derivable from the sum of the two potentials, 
4>a{p,z) = (l)ai{z) + (t>a2{p,z), Or the sum of the two 
stream functions, ipa{p, z) = ipaiip) + 4'a2{p,z), via 
Haz = -d(t)a{p,z)/dz = -~{l/p)d'ipa{p,z)/dp and Hap = 
—d(j)a{p,z)/dp= {1/ p)dipa{p,z)/dz. In the plane of the 
disk^ 



Haz{p,0) = 



0, p < a, 

TT V 

p> a, 



(A14) 



H. 



a -1 ft 

= — sm - 



(A15) 



and Hap{p,0+) — i?Q2p(p, 0+), which is given in Eqs. 
(IM) and (IXm . Since Hapip,0-) = -Hap{p,0+), the 
induced sheet-current density in the disk is {p < a) is 



Ka4,ip)=2Hap{p,0+) 



It is easily shown that 



4.HaP 



TT^ — p^ 



(A16) 



/ [Ka^{p)l2p\dp = Ha2z{0, 0) = -Ha, (A17) 

Jo 



as required by the Biot-Savart law. 



APPENDIX B: MAGNETIC FLUX <I> AT THE 
DISK'S CENTER 



Consider next the magnetic field and induced current 
density when a vortex carries magnetic flux $ = c^o = 
h/2e up through the center of the disk, or when the di- 
ameter 6 of a flux dome carrying magnetic flux $ = N(j)Q 
up through the disk is much smaller than the radius a 
of the disk. Neglecting magnetic structure in the dome 
on the length scale of b, A, or A, the z component of the 
magnetic field generated by the vortices is 



Hyz{p,0) 



6{p), p<a, 
<i>a 



(Bl) 
p> a. (B2) 



Both Hyzip,z) and the corresponding radial compo- 
nent Hyp{p,z) can be expressed in terms of the mag- 
netic potential 4>y{p,z) and stream function ijjy{p, z) via 
Hyz{p,z) = -d(t)y{p,z)/dz = -{1/ p)dipy{p,z)ldp and 
Hyp{p,z) = -dcpy{p, z)/dp ^ {l/p)dijy{p, z)/dz, where 



: 0. The solutions for z > are 

POO 

(t>v{p,z) = / fy{k)jQ{kp)e~^''dk 



■^v{.P,z) = 

Hyz{p,z) = 

Hyp{p,z) = 

Mk) = 

g(ka,0) = 

si(z) = 
In the plane z — 0, 

Mp. 0+) 



fy{k)Jo{kp)e-'''-kdk 



(B3) 

fy{k)Ji{kp)pe-^'dk, (B4) 
(B5) 
(B6) 
(B7) 

dp = si(fca) (B8) 
(B9) 



fy{k)Ji{kp)e~''''kdk, 
—g{ka,0) 



27r^o P-o 

aJojkp) 

a p^ p^ - a2 

sini 



t 



-dt. 



* -iP ^ 
cos — , p < a, 

a 



TT^POP 



0, p > a. 



^v{p, 0- 



Hvpip,0- 



$ TT 



tt-'Pq 



sin — , p> a, 
P 

( P 



„ „ I , , cos ^ - 



0, p > a, 



(BIO) 

(Bll) 

), P <a, 
(B12) 



and the sheet-current density in the disk (p < a) is given 
by Ky^{p) — 2Hyp{p,0+). The magnetic field hues gen- 
erated by the flux $ at the origin are illustrated in Fig. 
[l] which shows contours of the stream function ipv{p, z). 



APPENDIX C: SMALL FLUX DOMES 

For the case of a small flux dome of radius h much 
smaller than the radius a of the disk, it is possible to 
calculate the dome-generated magnetic fields Hsz{p,z) 
and Hsp{p, z) as well as the sheet-current density Ks^{p) 
generated by vortices in the dome. Consider the case of a 
disk of infinite radius, where the magnetic field generated 
by the vortices is 



Hsz{p,0) - Hd^/l-p^b^, p<b, (CI) 
= 0, p> 6. (C2) 

Both Hsz (p, z) and the corresponding radial compo- 
nent Hsp{p,z) can be expressed in terms of the mag- 
netic potential (/)s {p, z) and stream function ips {p, z) via 
Hsz{p,z) = -d(j)s{p,z)/dz = -{l/p)dips{p,z)/dp and 
Hsp(p,z) = -d(j)s{p,z)/dp = {l/p)dtps{p,z)/dz, where 
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W^(j)s = 0. The solutions for 2; > are 

/■OO 

(j^,{p,z) = / fs{k)Jo{kp)e-^'dk, (C3) 
Jo 

/•OO 

V«(p,2) = - / fs{k)Jiikp)pe-'''dk, (C4) 
Jo 



POO 

Hsz{p,z) = / f,{k)Jo{kp)e-^'kdk, 
Ja 



Hsp{p,z) = 

Jo 

fs{k) = Hdb^ 
In the plane z = 0, 



/ fs{k)Ji{kp)e-'''kdk, 
Jo 

sin kb — kb cos kb 



{kbf 



(C5) 
(C6) 
(C7) 



(C8) 



p> 6. 



V'a(P,0+) 



3 

3 



, P<b, 



H,p{p,0+) = Hd-^, p<b, 



4 6' 



(C9) 

, p>6,(C10) 
(Cll) 



p . -lb y/p"^ - 62 

-sin 

bp p 



P>b, 

and the sheet-current density is given by 

Ks4.{p) = 2Hsp{p,0+). 
It can be shown that 

/ [Ks^{p)/2p]dp = H,,{0, 0) = Hd 
Jo 

as required by the Biot-Savart law. 



(C12) 
(C13) 

(C14) 



APPENDIX D: MAGNETIC FLUX # UP 
THROUGH AN ANNULUS OF RADIUS pc 

When a closely spaced array of vortices carries mag- 
netic flux $ up through a narrow annulus of radius 
Pc < a, the magnetic field and induced current den- 
sity the vortices generate can be calculated as follows. 
Neglecting magnetic structure on the length scale of the 
intcrvortcx spacing, A, or A, the z component of the mag- 
netic field generated by the vortices in the plane z = 
is 



Hc.{p,0) 



2niiop 



Hp- Pc), P<a, 



7rVo(p^ - pDVp'^ - «^ 



(Dl) 
p > a.(D2) 



Both Hcz{p,z) and the corresponding radial compo- 
nent Hcp{p, z) can be expressed in terms of the mag- 
netic potential (j)c{p, z) and stream function ipcip, z) via 
Hcz{p,z) = -d(j)c{p,z)/dz = -{1/ p)d'4}c{p,z)ldp axid 
Hcp{p,z) = -d4>c{p,z)ldp = {l/p)dtpc{p,z)/dz, where 
V^0c = 0. The solutions for 2; > are 



Mp^z) 
i)c{p,z) 

Hcz{P,z) 

Hcp{p,z) 
fc{k) 
g{ka,pc/a) 



/ fc{k)Jo{kp)e-^'dk, (D3) 
Jo 

f-OO 

- / fc{k)Ji{kp)pe-'''dk, (D4) 
Jo 

POO 

/ fc{k)Jo{kp)e-^'kdk, (D5) 
Jo 

PQO 

/ fc{k)Ji{kp)e-'''kdk, (D6) 
Jo 



2ttpo 

--f 

J a 



Jo{kpc 



TT^PO 



g{ka,pc/a)lD7) 



P\Ja^ - plJojkp) 
(p2 - pD^fP -a? 



dp. (D8) 



In the plane 2 = 0, 

(Ac(p,0+) = 

V'c(p,0+) = 



■^{a,Pc,p), P<a, 



2npo 
0, p> a, 

0, p < Pc, 

$ TT 

2 77, Pc<P<a, 



(D9) 



TT^PO 



sm 



P^ -pi 



, p > a,(D10) 



Hcp{p,0+) = 



2ttpo 
= 0, p > a 



G{a,pc,p), p<a, 



(Dll) 



where, defining p< (/!?>) to be the smaller (larger) of pi 
and p2, 

T{a,pc,p) = F{pc,p) 





2 r ^a^^plF{p\p)p' 
T^Ja {p'^ ~ pD^p'^ -a^ 


V,(D12) 


F{Pl,P2) 


POO 

= 1 Jo{kpi)Jo{kp2)dk 
Jo 


(D13) 




= —K{p^/p^), 
7rp> 


(D14) 


F{p,0) 


= F{0,p) = l/p, 


(D15) 


Ha,p',0) 


2 -1 p' 
= — -cos — , 
wp' a 


(D16) 


T{a,p',a) 


= 0, 


(D17) 
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G{pi,P2) = 



G(0,p) 
G{p,Q) 



-dT{a,pc,p)/dp 
G{pc,p) 
2 

TT 



V^^~7cG{p\p)p' 



-dF{pi,p2)/dp2 



dp', 



Jo{kpi)Ji{kp2)kdk 



(-) 



K{k) E{k) 



-P1 + P2 Pi 
/ 2 X E{p,/p2) 
\Trp2J P2{1- pI/pD' 



(-) 

Pi + P2 ' 

0, 



JS:(P2/pi) 



- P2 
Pi < P2, 

E{p2/pi) 



(D18) 
(D19) 
(D20) 

(D21) 

(D22) 



Pi 



pi{i-pI/pI)- 

pi > p2, (D23) 
(D24) 

(D25) 
(D26) 



Jo{x) and Ji{x) are Bessel functions, and K{k) and E{k) 
are complete elliptic integrals of the first and second kind 
of modulus K. The sheet-current density in the disk {p < 
a) is given by Kc,f,{p) = 2Hcp{p,0+). It can be shown 
with the help of Eqs. (|D5]) and jHl) that 



[K,^{p)/2p]dp = H,,{0,0)=Q, 



(D27) 



as required by the Biot-Savart law. A plot of the mag- 
netic field lines for this case would be similar to that 
shown in Fig. [TJ except that the field lines would em- 
anate from p = pc rather than p = 0. 



^ J.-P. Girard, E. Paumier, and A. Fortini, J. Phys. (Paris) 
32, 895 (1971). 

^ R. P. Huebener, R. T. Kampwirth, and J. R. Clem, J. Low 

Temp. Phys. 6, 275 (1972). 
^ J. R. Clem, R. P. Huebener, and D. E. Callus, J. Low 

Temp. Phys. 12, 449 (1973). 

* J. Provost, E. Paumier, and A. Fortini, J. Phys. F 4, 439 
(1974). 

^ H. Castro, B. Dutoit, A. Jacquier, M. Baharami, and L. 

Rinderer, Phys. Rev. B 59, 596 (1999). 
^ M. V. Indenbom, H. Kronmiiller, T. W. Li, P. H. Kes and 

A. A. Menovsky, Physica C 222, 203 (1994). 

Th. Schuster, M. V. Indenbom, H. Kuhn, E. H. Brandt, 

and M. Konczykowski, Phys. Rev. Lett. 73, 1424 (1994). 

* E. Zeldov, A. I. Larkin, V. B. Ceshkenbein, M. Kon- 
czykowski, D. Majer, B. Khaykovich, V. M. Vinokur, and 

H. Shtrikman, Phys. Rev. Lett. 73, 1428 (1994). 

A. Fortini and E. Paumier, Phys. Rev. B 14, 55 (1976). 
A. Fortini, A. Haire, and E. Paumier, Phys. Rev. B 21, 
5065 (1980). 

" M. Benkraouda and J. R. Clem, Phys. Rev. B 53, 5716 
(1996). 

^2 T. B. Doyle, R. Labusch, and R. A. Doyle, Physica C 290, 
148 (1997). 

I. L. Maksimov and A. A. Elistratov, Appl. Phys. Lett. 72, 
1650 (1998). 

" M. Benkraouda and J. R. Clem, Phys. Rev. B 58, 15103 
(1998). 

Y. Mawatari and J. R. Clem, Phys. Rev. Lett. 86, 2870 
(2001). 

N. V. Zhelezina and G. M. Maksimova, Tech. Phys. Lett. 
28, 618 (2002). 

A. A. Babaei Brojeny, Y. Mawatari, M. Benkraouda, and 
J. R. Clem, Supercond. Sci. Technol. 15, 1454 (2002). 
V. A. Schweigert and F. M. Peeters, Phys. Rev. Lett. 83, 
2409 (1999). 



P. Singha Deo, V. A. Schweigert, and F. M. Peeters, Phys. 

Rev. B 59, 6039 (1999). 
2° E. H. Brandt, Phys. Rev. B 60, 11939 (1999). 
2^ E. H. Brandt, Physica C 332, 99 (2000). 
22 E. H. Brandt, Low Temp. Phys. 27, 723 (2001). 

A. A. Babaei Brojeny and J. R. Clem, Phys. Rev. B 68, 

174514 (2003). 

'^^ G. R. Berdiyorov, L. R. E. Cabral, and F. M. Peeters, J. 

Math. Phys. 46, 095105 (2005). 
2^ J. Pearl, Appl. Phys. Lett. 5, 65 (1964). 
2® P. N. Mikheenko and Yu. E. Kuzovlev, Physica C 204, 229 

(1993). 

J. R. Clem and A. Sanchez, Phys. Rev. B 50, 9355 (1994). 
M. Yu. Kupriyanov and K. K. Likharev, Fiz. Tverd. Tela 
16, 2829 (1974) [Sov. Phys. Solid State 16, 1835 (1975)]. 
N. Morozov, E. Zeldov, M. Konczykowski, and R.A. Doyle, 
Physica C 291, 113 (1997). 

A. V. Kuznetsov, D. V. Eremenko, and V. N. Trofimov, 
Phys. Rev. B 56, 9064 (1997). 

A. V. Kuznetsov, D. V. Eremenko, and V. N. Trofimov, 
Phys. Rev. B 57, 5412 (1998). 

Y. M. Wang, A. Zettl, S. Ooi, and T. Tamegai, Phys. Rev. 
B 65, 184506 (2002). 

L. G. Aslamazov and S. V. Lempicki, Zh. Eksp. Teor. Phys. 
84, 2216 (1983). 
^'^ E. Zeldov, J. R. Clem, M. McElfresh, and M. Darwin, 
Phys. Rev. B 49, 9802 (1994). 

D. Yu. Vodolazov, I. L. Maksimov, and E. H. Brandt, Eu- 
rophys. Lett. 48, 313 (1999). 

R. B. Doyle, R. Labusch, and R. A. Doyle, Physica C 332, 
365 (2000). 

A. A. Elistratov, D. Yu. Vodolazov, I. L. Maksimov, and 
J. R. Clem, Phys. Rev. B 66, 220506(R) (2002); erratum 
Phys. Rev. B 67, 099901(E) (2003). 
^® Y. Mawatari and J. R. Clem, Phys. Rev. B 68, 024505 
(2003). 



10 



J. D. Jackson Classical Electrodynamics (Wiley, New York, ^° H. Lamb, Hydrodynamics (Dover, New York, 1945), p. 138. 
1962), p. 146. 



